clc
clear
for j=1:3
mpc=loadcase(casomasgrandebucle4)
rundcopf(casomasgrandebucle4) % resuelve el opf.
define_constants;
resultado=rundcopf('casomasgrandebucle4');
Mult_Lagrange(:,j)=resultado.bus(:,LAM_P);
end
display(Mult_Lagrange)